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We present an analysis of Brownian force noise from residual gas damping of reference test masses 
as a fundamental sensitivity limit in small force experiments. The resulting acceleration noise 
increases significantly when the distance of the test mass to the surrounding experimental apparatus 
is smaller than the dimension of the test mass itself. For the Advanced LIGO interferometric 
gravitational wave observatory, where the relevant test mass is a suspended 340 mm diameter 
cylindrical end mirror, the force noise power is increased by roughly a factor 40 by the presence 
of a similarly shaped reaction mass at a nominal separation of 5 mm. The force noise, of order 
20 fN/Hz 1 / 2 for 2 x 10 -6 Pa of residual H2 gas, rivals quantum optical fluctuations as the dominant 
noise source between 10 and 30 Hz. We present here a numerical and analytical analysis for the gas 
damping force noise for Advanced LIGO, backed up by experimental evidence from several recent 
measurements. Finally, we discuss the impact of residual gas damping on the gravitational wave 
sensitivity and possible mitigation strategies. 

PACS numbers: 05.40.-a, 07.10.Pz, 07.30.-t, 95.55.Ym 



I. INTRODUCTION 

The Brownian motion of a macroscopic test mass is 
a fundamental limit in many small force measurements, 
which ideally require a test mass (TM) to be free of stray 
forces, and thus a reference of purely inertial motion, ex- 
cept for any known and calibrated suspension forces. The 
power spectrum of the Brownian force noise acting on the 
TM is related to any source of mechanical dissipation, 
given by the mechanical impedance Z (oj) , through the 
fluctuation-dissipation theorem 



suiting force noise, over that observed for the same TM in 
an infinite volume filled only with gas. This phenomenon, 
referred to here as proximity-enhanced gas damping, has 
also been observed and studied for MEMS oscillators, 
under the name of squeeze film damping [7|, [8| . 

While there is a continuum of behavior between the 
free damping (infinite gas volume) limit and proximity 
damping, we will divide the total molecular impact force 
noise into these contributions, 



S F (w) = Sf + AS F (w) 



(2) 



S F (w) = 4k B TRe [Z [oj)\ = 4k B TRe 



Fjto) 



(1) 



Gas damping in the molecular flow regime is known to 
produce viscosity proportional to the residual gas pres- 
sure, Z (lo) = j3 ~ pjjHJI, with f3 referred to here as 
the gas damping coefficient. In the limit to which the 
collisions can be treated as independent impulses, the 
resulting Brownian force noise has a frequency indepen- 
dent spectrum. Residual gas force noise has recently been 
reconsidered in the context ofgravitational experiments 
with geodesic-reference TM 00]. In these and other 
experiments, the proximity of surrounding surfaces can 
significantly increase the gas damping, and thus the re- 



Sf = ^ksTp 00 is the noise for the TM in an infinite 
gas volume, while the excess ASV depends on the TM 
proximity to surrounding surfaces. 

Though long unnoticed in experimental gravitation, 
proximity-enhanced gas damping has a straightforward 
physical explanation: motion of a TM, with character- 
istic dimension s, in the vicinity of another surface, at 
some small distance d, creates a transient squeezing of 
any residual gas in the gap between the two bodies. With 
the subsequent molecular flow out of the gap, there is an 
associated pressure drop, proportional to the molecular 
current, and thus a velocity dependent force. In an alter- 
nate picture, the Brownian force noise increases because 
of the correlation between repeated impacts of a molecule 
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on the TM as it stochastically moves along the gap. As 
the same random walk statistics govern the molecular 
flow down the channel, these two pictures - of dissipative 
flow and correlated collision impulses - are equivalent, as 
Eqn. [T] requires. 

The grouping of repeated impulses into "macro- 
collisions" establishes a natural interaction (or correla- 
tion) time r, the typical time needed for random walk 
molecular diffusion out of the interbody gap. This cre- 
ates a high frequency cutoff, to t _1 , above which the 
force noise decreases. In terms of the macroscopic flow 
impedance picture, TM motion for lot ^> 1 is too fast to 
allow molecular flow, with the real, dissipative, part of 
Z (uj) decreasing to leave only the imaginary impedance 
of gas compression. We will return to these simple argu- 
ments later in the article for an approximate analytical 
model that is useful in interpreting the results of our 
study, for ASf and r in terms of the aspect ratio s/d. 

Recent experimental studies [g, |6j place the physics of 
proximity-enhanced gas damping on a quantitative foot- 
ing. These studies both measure the gas damping on the 
free oscillation of torsion pendulums, focusing on the low 
frequency range lot <C 1 where ft is independent of fre- 
quency. The torsion pendulum inertial elements in these 
studies were a cubic TM inside a cubic enclosure and 
a rectangular plate in close proximity of a second paral- 
lel plate. The TM size-to-gap aspect ratios studied range 
from roughly 4 to 400. Both experiments observed a clear 
proximity excess in the gas damping, well beyond calcu- 
lated values of /3°° . These experiments will be discussed 
in Sec. [TV] 

Gas damping has been analyzed as a force noise source 
for terrestrial gravitational wave interferometers [3| but 
should be reconsidered for the next generation experi- 
ments, both for their stringent acceleration noise require- 
ments and for the close proximity of TM to the surround- 
ing apparatus. For Advanced LIGO [H, the envi- 
sioned end mirror TM are right cylinders of fused silica 
with radius R = 17 cm and length h — 20 cm (mass 
M w 40 kg). In an infinite gas volume, at T = 293 K 
and pressure p = 2x 10~ 6 Pa of residual H2 gas (mo = 
2 amu), the molecular impact force noise, integrated over 
the entire cylindrical TM surface, would be [Jj 



■V = = !>ir2SKnU Li T) 1/2 R 2 (l + ^ + J 



3.2fN/Hz 



1/2 



P 



2 x 10" 6 Pa 



1/2' 



(3) 



This free damping compares favorably with the Ad- 
vanced LIGO target, which requires total force noise be- 
low 40 fN/Hz 1 / 2 at frequencies near 30 Hz to reach a 
strain sensitivity of 10~ 23 /Hz 1 ^ 2 (equivalent to a TM ac- 
celeration noise of roug hly 1 fm/s 2 /Hz 1/2 ). However, the 
original Advanced LIGO design includes cylindrical "re- 
action masses," with the same radius and 13 cm length, 
facing the TM at a nominal separation d = 5 mm (see 



2R K 



He 



RM 



TM 

h 



2R 




FIG. 1: Cartoon illustration of the Advanced LIGO exper- 
imental configuration studied and simulated here, seen per- 
pendicular to the symmetry axes of the cylindrical test (TM) 
and reaction (RM) masses, with radii R = 17 cm and heights 
h = 20 cm and Hrm = 13 cm. The gap d, nominally 5 mm for 
Advanced LIGO, is varied in the simulations. The simulated 
experimental enclosure, not shown to scale here, is a coaxial 
cylinder of radius Rs = 10R and length Hs = Uh + fiRM +d. 



Fig. [T]) . The close reaction masses are needed for all four 
cavity mirror TM, for electrostatic force actuation and 
for compensating optical thermal lensing from TM laser 
heating. The relevant aspect ratio of 2R/d « 70 is well 
into the regime of proximity-enhanced gas damping, and 
thus the problem merits a thorough understanding. 

In this article, we calculate the force noise and diffusion 
time using a numerical simulation of a molecule exchang- 
ing momentum with the TM mirrors in the Advanced 
LIGO geometry, as a function of the gap size (Section 
ITT|) . In Section Hn] these results are compared with an ap- 
proximate analytical model, which successfully explains 
the functional dependence of the excess force noise on gap 
size in terms of a simple random walk picture. In Section 
IIVI the physics and simulation techniques of Sees. ILTland 
IIIII are compared with experimental evidence from simi- 
lar systems. Finally, in Section [V] we address the impact 
of gas damping force noise on the sensitivity of advanced 
gravitational wave detectors and discuss possible mitiga- 
tion strategies. While the article focuses on the Advanced 
LIGO geometry, we will discuss design consideration for 
other advanced gravitational wave observatories in the 
conclusion. 

As in previously cited experimental and numerical 
studies 01,1 , we consider completely inelastic molecu- 
lar collisions with the TM, with subsequent immediate 
random reemission with a cosine angular distribution. 
The inelastic, diffuse scattering hypothesis is verified ex- 
perimentally for spinning rotor pressure gauges and is the 
basis for their use as absolute pressure calibration stan- 
dards [III . This also guarantees time reversal invari- 
ance and Maxwcll-Boltzmann statistics in the gas phase 
(see Ref . [l3| ) . The immediate reemission approximation 
is valid in the limit that the effective sticking time, r st , is 
negligible compared to that needed to cross the gap be- 
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tween TM, This is justified for H2 and other possibie 
gas species on the TM SiC>2 surface at room temperature, 
based on the relevant adsorption potentials, in Appendix 
|A1 Finally, we treat the average gas pressure as uniform 
in the system, which requires that the gap pressure is not 
dominated by local outgassing from the test or reaction 
mass, which would cause a stable increase in the relevant 
local pressure. We return to these last two points in the 
conclusions. 



II. NUMERICAL STUDY OF BROWNIAN 
FORCE NOISE 

We have performed Monte Carlo simulations to ana- 
lyze the gas damping noise for the Advanced LIGO geom- 
etry. The simulations (following those performed in Q) 
trace the trajectory of a gas molecule as it moves inside a 
large volume containing the cylindrical test mass and re- 
action mass. The molecular species is chosen to be H2 at 
room temperature, the dominant residual gas expected 
in the Advanced LIGO vacuum chamber, which sets the 
scale of momentum exchange upon impact and the time 
scale between collisions through the characteristic ther- 
mal velocity vt = yk^Tjm ~ 1100 m/s. The volume 
enclosing the two cylinders, ideally infinite, is chosen to 
be a cylinder ten times larger, in both radius and length, 
than the volume envelope defined by the two cylinders 
and the intervening gap. 

The simulation starts by random selection of a molecu- 
lar position, from a uniform distribution in the available 
volume - including the gap but also the much larger, 
and thus more probable, space surrounding the cylin- 
ders - and velocity, from a Maxwell-Boltzmann distribu- 
tion. The impacts with the walls are completely inelastic, 
with subsequent immediate random thermal reemission 
following a cosine angular distribution, independent of 
the incoming velocity, with the resulting gas distribu- 
tion obeying Maxwell-Boltzmann statistics (the reemis- 
sion distribution is given in Eqn. ICll as in Ref. Q). At 
each collision with the TM -regardless of whether the 
molecule strikes on the TM surface facing the gap or on 
the opposing side or outer wall - the time and exchanged 
momentum vector are recorded. Simulations are run, for 
a single molecule, for a time To (in the range of 10~ 6 
- 10 s, chosen to abundantly cover the full range of the 
observed values of t) and then repeated for many trials. 
Finally, the simulations have been performed as a func- 
tion of gap d between the two cylinders, in the range from 
100 (an to 1 m. 

One way to probe the simulation is to observe the 
statistical fluctuations of the total momentum Aq ex- 
changed along the critical cylindrical symmetry axis in 
repetitions of the simulations of duration To. This allows 
simple extraction of the parameters describing the ampli- 
tude and spectral shape of the force noise. For a process 
with white force noise up to some high frequency rolloff 
associated with a time constant r, with Sf = , , , F ° ^ , 
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FIG. 2: Simulation data for the scatter in the total momen- 
tum Aq transferred to the TM in simulation time To, for 
several values of the intercylinder gap d, each fit to the model 
in Eqn. [6] The range of gaps studied represents aspect ratios 
(2R/ D) from nearly 7 to 340. The dashed gray line repre- 
sents the prediction in the limit of infinite gap, S^, which is 
recovered for very short simulation times. Data are scaled for 
a reference pressure p — 2 xlCP 6 Pa. 

one can calculate (see Appendix iBj) 

«*,)■> -^(l- 1 , (4) 

with ((Aqf) = in the limit T » r. 

In our simulations, we directly calculate the average 
scatter in Aq from the total momentum exchange in the 
different trials, typically calculating ((Aq) 2 ) from groups 
of 10 3 (for the longer simulations) to 10 6 (shorter To) 
single molecule experiments and then estimate the simu- 
lation uncertainty between 100 such groups. In a physi- 
cal experiment with many non-interacting gas molecules, 
((Aq) ), and thus also the force noise power, is propor- 
tional to the number of molecules in the system, allowing 
us to scale our single molecule simulations to a given pres- 
sure (we choose 2 x 10~ 6 Pa as a reference pressure to 
put a relevant experimental scale to the data presented) . 

Figure [2] shows simulation data for '- as a func- 

tion of measurement time To for several vafues of the gap 
d. The data show a clear increase in force fluctuations 
with decreasing gap. The time dependence shows satu- 
ration for large To, as predicted by Eqn. [4] and which we 
can associate with r, but also a saturation to a nonzero 
value at short To, for which this simple model with a sin- 
gle high frequency rolloff does not account. This low- To 
saturation has, however, a simple physical explanation: 
for To below roughly d/vx - the time necessary for a 
molecule to cross the intercylinder gap, see Appendix [C] 
- a molecule will hit the TM at most once during the sim- 
ulation. In this limit the simulation becomes blind to any 
correlation effects related to proximity, as for a cylindri- 
cal TM in an infinite volume. In fact, the low- To satura- 
tion values shown in Fig. [5] converge, for different values 
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of d, to the infinite volume value, Sp?. There would also 
be a natural time scale and associated high frequency 
rolloff relevant to the simulation in the single collision 
(or infinite volume) limit, but it corresponds to the sin- 
gle collision duration, which the simulation takes to be 
infinitesimal. These curves thus both confirm the infinite 
volume force noise prediction for short simulations and, 
in the longer simulations, indicate an increasing excess 
force noise for decreasing d. 

We note that no excess force fluctuations are observed 
for the components orthogonal to the cylinder symme- 
try axis, for which the infinite volume limit of Ref. Q 
is verified. The correlation in the force from successive 
impacts is only non-zero for the component normal to 
the surface, for which the direction of the impulse is al- 
ways toward the test mass center, whereas the sign of the 
shear component of momentum exchange along an axis 
parallel to the surface is random. From a macroscopic 
damping standpoint, TM motion parallel to the reaction 
mass does not squeeze the gas between the two bodies, 
and thus creates no dissipative gas flow in the narrow 
gap. 

With these results in mind, we modify the force noise 
model to 

S F (u,) = Sf + - Fo (5) 
1 + (ur) 

where ASf represents the excess low frequency force 
noise associated with the proximity effect. For each d, 
the curves in Fig. [2] are thus fitted to a model modified 
from Eqn. [5] to include the infinite volume limit, 

with the total low frequency force noise Sp obtained in 
the limit To 3> r, 

T -»-oo 1q 

This allows extraction of the free fit parameters corre- 
sponding to the low frequency force noise Sf and dif- 
fusion time r, which are plotted as a function of d in, 
respectively, Figs. [3] and 2] [28[ . 

The Monte Carlo molecular dynamics simulation can 
also be probed by sampling the momentum exchange 
with the test mass in repeated time intervals to create 
a force time series, which can then be analyzed for the 
frequency dependent noise spectrum. Force noise spectra 
have been calculated using a 100 fis sampling time - and 
thus well shorter than the relevant values of r - for a sim- 
ulation of 10 seconds, which has been analyzed in over- 
lapping 0.4 s windows, with averaging over roughly 25 
statistically independent spectra. These data, displayed 
for 5 and 20 mm gaps in Fig. [5J are in agreement with 
the model in Eqn. [5j which has been plotted on top of 
the noise spectrum data using the values of ASf an d t 
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FIG. 3: Simulation data (blue dots) for the force noise lin- 
ear spectral density as a function of the gap, d, for p — 2 
xlO -6 Pa, with a fit to the simple analytical model (Eqn. 
I20[) . The dashed gray line is the theoretical infinite gap limit. 
Points at d of 5 and 20 mm are circled to indicate, respec- 
tively, the geometries of the baseline and alternative designs 
discussed in the conclusion. Points shown as black "x" are 
simulation data for the molecular escape time, converted into 
force noise using the model in Sec. IHI Al 
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FIG. 4: Value of the characteristic time constant r extracted 
from the simulation force fluctuations (blue dots), as a func- 
tion of gap d, using the model in Eqn. [6] The dashed pink 
curve represents the raw result of the approximate analytic 
model (Eqn. [2} in Sec lIII Bl while the solid blue curve repre- 
sents a two parameter fit to the same model, discussed in the 
text. Again, circled points at 5 and 20 mm indicate the base- 
line and alternate geometries. Points marked with a black 
"x" are simulation values for the average molecular escape 
time from the intercy Under gap, which correspond approxi- 
mately to the values of r extracted from the force fluctua- 
tions. 



obtained from {(Aq) ) as described above and presented 
in Figs. CU-H 

A variation on this simulation can be used to extract 
the average escape time for a molecule inside the inter- 
cylinder gap. Here, the initial position of the molecule 
is distributed uniformly inside the cylindrical space be- 
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Frequency (Hz) 

FIG. 5: Simulation force noise linear spectral density as a 
function of frequency for 5 and 20 mm gaps, using p = 2 x 
10 -6 Pa. Also plotted, with solid black lines, is the model 
from Eqn. [5l using the parameters extracted from the data in 
Fig. [2] The dashed gray line corresponds to the force noise in 
the infinite volume limit. 



tween TM and RM, with the simulation run until the 
molecule exits this gap. This has been performed for 
7 values of d, each time averaging over 10000 trial 
molecules. These data are plotted with the values of 
t extracted from the force fluctuations in Fig. 01 The 
agreement, at the 10-15% level, with the time constant 
relevant to the force noise spectrum rolloff, confirms the 
intuitive interpretation of diffusion time discussed in the 
introduction. Additionally, as will be discussed in the fol- 
lowing section, the diffusion time t (d) is linked directly, 
with a simplified analytical model, to an estimate of the 
resulting force noise. This projection for the white force 
noise level is shown here in Fig. [3] 

Before discussing the quantitative impact of the ob- 
served noise increase on the gravitational wave detection 
sensitivity, we turn first to a discussion of the physics of 
excess gas damping in the context of simplified analytical 
models and then to the experimental evidence underpin- 
ning the simulation and predictions. 



III. COMPARISON WITH SIMPLIFIED 
ANALYTICAL MODELS 

A combination of simple analytical models is useful 
for illuminating the physics of proximity-enhanced gas 
damping, but also for arriving at a closed form approx- 
imate expression for the excess force noise as a function 
of the gap d. 



A. Macroscopic damping for a piston weakly 
connected to an external gas reservoir 

A simple macroscopic model for proximity-enhanced 
jas damping, following the discussion in the introduction 



and Refs. Q, comes from approximating the inter- 
cylinder gap as an isothermal volume at uniform pressure 
connected to the external gas reservoir at pressure p via a 
molecular flow impedance, Zm, with units of [s/m 3 ]. In 
the presence of a displacement x (t) of the TM along its 
symmetry axis, the volume is V — irR 2 (d + x), and we 
can express the pressure as p + Sp (t) . We can then write 
the standard ideal gas and molecular flow equations for 
the instantaneous number of molecules Nq in the gap, 

k B TN = (p + Sp) ttR 2 (d + x) 

k BT—r- = -— (8) 
at Zm 

Considering small displacements i < d, these expres- 
sions can be combined into a linearized equation for Sp, 

dSp Sp _ _P<^_ (( .\ 
dt t d dt 1 ' 

where we have reexpressed the molecular impedance in 
terms of the effective pressure equilibration time r = 
ZmV = ZMirR 2 d, analogous to an RC discharge time. 

To calculate the mechanical impedance, or relationship 
of force to velocity v, Eqn. [§] can be evaluated in the 
frequency domain for harmonic TM motion, using v = 
iux. The pressure increase Sp is multiplied by surface 
area to obtain force, 

2 pr 1-iuiT 



F = -ttR 



2 pT ILOX 

d 1 + iuj 



= -v X ttR 



(10) 



d 1+ {lot) 2 

Equation [10] describes a mechanical impedance that is 
primarily real and thus dissipative at lower frequencies. 
At higher frequencies, the impedance becomes imaginary, 
and thus a lossless spring as gas compression occurs with- 
out sufficient time for gas flow out of the gap. For the 
force noise, the fluctuations-dissipation relation in Eqn. 
[T] can be evaluated with the real part of Eqn. [TU] to yield 

1 



AS F = 4:k B T x ttR 



2P T 



(11) 



d 1 + {lot) 2 

This low-pass form of the predicted force noise is indeed 
confirmed by the simulations, in Figs. [2] and [5] 

For the interaction-free regime studied here, the pres- 
sure equilibration time r is equivalent to the average time 
that a molecule needs to escape from the gap into the 
surrounding gas reservoir. This escape time is plotted 
in Fig. U and, converted into a low frequency force noise 
level using Eqn. QTJ Fig. [3] This macroscopic picture is 
approximate in that it assumes the pressure inside the 
gap to be uniform, while physically there must be a ra- 
dial pressure gradient in order to drive gas flow out of 
the gap upon squeezing. Nonetheless, there is a remark- 
able agreement for the cylindrical geometry between the 
average escape time and the simulation force fluctuation 
values for force noise and time constant r. 

To complete a predictive analytical model for proxim- 
ity damping, we require an independent estimate for the 
time constant r (d) , which contains the physics of the 
molecular diffusion process and is addressed in the fol- 
lowing subsection. 
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B. Random walk calculation of escape time r 

The escape time r can be estimated by simple random 
walk arguments for a molecule in the intercylinder gap. 
This model is relevant in the limit d <C R, with a typical 
molecule making many collisions before escaping the gap. 

The number of collisions N that a particle makes be- 
fore diffusing out of the intercylinder gap can be esti- 
mated from the random walk mean square distance trav- 
eled along a direction perpendicular to the cylinder axis, 
or N(r 2 ), where (r 2 ) is the mean square lateral displace- 
ment that a particle makes in a single flight between op- 
posing faces of the two cylinders. Thus, covering a dis- 
tance R to diffuse out of the gap will require N ps t^jt. 
The characteristic time r that a molecule takes to ran- 
dom walk out of the gap will thus be r ps N(t), where 
(i) is the average time of flight to cross the gap between 
opposing faces of the two cylinders. 

Estimates of (r 2 ) and (i) can be calculated analytically 
for a particle emitted, with cosine-law distribution, from 
the center of one of the cylinder faces into the gap of 
height d and radius R (see Appendix [Cjl : 



(r 2 ) 



d 2 In 



V 2 / vt 



We can use these values to estimate N, 

R 2 



N 



R 2 ^ 

( r2 > ~ d 2 In 



(12) 



(13) 



and then the molecular escape time r, 
. 1/2 E 2 



dvT In 



(14) 



Inserting the result of Eqn. [14] into Eqn . ITT1 and recall- 
ing the thermal velocity vt = \/k B T /m, we obtain the 
low frequency force noise 



AS Fo ^p(8irm k B T) 1/2 TrR 2 - 



R 2 



d 2 In 



1 



(15) 



The force noise thus has a prefactor - proportional to p 
and surface area irR 2 - similar to the infinite volume re- 
sult (Eqn. [3]), multiplied by a dimensionless excess factor 
that depends on the aspect ratio via R/d. 

This model for r and Sf is also approximate. The true 
mean square displacement (r 2 ) depends on the position 
on the cylinder face from which the particle is emitted, 
as will the expected full escape time. Additionally, the 
model treats the random walk as a Gaussian process, 
with the mean square step size as the single parameter 



characterizing the statistics. For R/d — > oo this approx- 
imation is catastrophic, as the variance of the step di- 
verges (see App. O and Ref. Q). However, for the 
scope of estimating the diffusion time out of a finite ra- 
dius cylindrical volume, the calculation of the step size 
variance can be truncated to consider steps no larger than 
R. As seen in earlier studies Q and as will be shown in 
the next sections, this is a good approximation to the 
physics seen in the simulations. 



C. Equivalent shot noise model 

Independently from the macroscopic piston damping 
model in Sec. IIII A[ the calculation of r in Sec. IIIIBI can 
be used for direct estimate of the force noise using a shot 
noise model. Here, the "shot" of momentum corresponds 
not to the momentum exchange of a single molecular im- 
pact, but rather to the total momentum exchanged by 
a single molecule over N impacts with the TM before 
escaping from the intercylinder gap. 

The force shot noise can be expressed 



AS Fo = 2 (Sq) 2 A 



(16) 



The characteristic impulse Sq is the typical thermal mo- 
mentum exchange of a single collision (2moWr) multiplied 
by the number N of collisions that a molecule makes with 
the TM on a typical pass through the the gap, 



Sq p» 2NttiqVt 



(17) 



The characteristic rate of such summed impulses, A, is 
the rate at which molecules enter, and exit, the inter- 
cylinder gap. This is the expected number of molecules 
between the cylinders divided by the diffusion time r, 



A 



' x nR 2 d x - 



k B T 



(18) 



Finally, insertion of these values into Eqn. [16] gives the 
force shot noise, 



2 (4N 2 m k B T) x 



P ,R 2 d X - 

T 

R 2 



k B T 



—p (SnmQksT) 1 ^ 2 nR 2 

77 ^ d 2 In 



1 



{19) 



This shot noise prediction differs from that of Eqn. Q~5] 
by the numerical factor ^ and thus, at the approximate 
level of this calculation, can be considered equivalent to 
Eqn.rj5l 



D. Comparison with simulation results 

Using Eqn. [15] as a model for the excess gas damping- 
force noise, a fit of the simulation data for Sf to the 
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model 



IV. SUMMARY OF EXPERIMENTAL 
EVIDENCE 



Bp(8TTmok B T) 1/2 TTR 2 



S Fo (d) « ASf + 
R 2 



d 2 ln !+(#) 



I) (20) 



is shown in Fig. [3] Here, S|? is denned as in Eqn. |3] and 
the subtraction of unity in the second term is inserted 
to give a null contribution for d — > +oo, attributing all 
residual damping to the infinite volume term. We obtain 
A = 1.038 ± 0.004 and B = 0.78 ± 0.005. The value 
of the excess damping factor B is of order unity, indi- 
cating that that this simplified model indeed gives the 
right magnitude of the excess damping, though the ob- 
served agreement at the 25% level is perhaps better than 
expected given the approximations of the model. It is 
likely that the finite size of the simulated system is re- 
sponsible for the several percent excess observed for the 
infinite gap limit (coefficient ^4). 

While the fit of this approximate model is not par- 
ticularly good by statistical standards, with \ 2 ~ 12, 
this model nonetheless gives a useful analytical formula 
for Sf as a function of gap, with 30% accuracy in noise 
power spectral density across four decades of d. We note 
that fitting to a model that neglects the logarithmic con- 
tribution - thus considering the excess damping term as 
ASp oc d~ 2 - deviates from the simulation data by more 
than a factor 2 across the same range. Similarly unsat- 
isfactory results are obtained with other power law fits, 
which can work well over a given single decade in d but 
not wider ranges, with the logarithmic slope observed to 
change, steepening towards -2 with decreasing gap. This 
confirms the basic physics given in Sec. IIII Bl 

The values of r extracted from Eqn. [5] are shown in 
Fig. |U These are compared with values associated with 
the simulation mean escape time from the gap. Shown as 
a dashed curve is a simple, unsealed plot of the approx- 
imate analytical estimate of r from Eqn. [14] 2j|, which 
is within a factor 2 of the simulation results across the 
three decades of d for which we extract t. Also shown, 
with a solid curve, is a two parameter fit to the data, us- 
ing Eqn. [14] but allowing multiplicative prefactors both 
to the overall expression for r and as a scale factor for d 
in the logarithmic term. This gives a statistically good 
fit, with the relevant scale factor effectively lengthening 
d in the logarithmic term, by 3.2 (equivalently, shorten- 
ing R). This is likely related to the oversimplification of 
assuming molecules starting from the geometric center 
- and not a range of positions in the intercylinder gap 

in our calculations of N, (t), and (r 2 ). We note that 
similar scaling of d in the logarithmic term for the force 
noise Sf 01 discussed in the preceding paragraph, has no 
relevant impact on the quality of the fit in Fig. |3] 



Two recent experimental efforts provide a verification 
of the general physics of proximity-enhanced gas damp- 
ing and of the simulation techniques employed in this 
study. Aside from the specifics of the geometries studied, 
the simulations used to analyze these experiments differ 
from those described in Sec. |TT] only in that they esti- 
mate a rotational damping coefficient, /3 rot — , using 
the fluctuations in the average torque, rather than force. 
Both experimental studies calculate the gas contribution 
to [3 ro t from measured torsion pendulum free-oscillation 
decay, characterized by energy decay time constant r e or, 
equivalently, the quality factor Q, 



1 



firot 



(21) 



where /o is the pendulum resonant frequency and / is its 
moment of inertia. 

Figure [5] displays the data for the measured ro- 
tational gas damping coefficient f3 rot from the Uni- 
versity of Trento (UTN) torsion pendulum ringdown 
measurements [5j with two pendulum geometries featur- 
ing cubic TM (46 mm side). Proximity-enhanced gas 
damping is studied in the geometry relevant to the LISA 
gravitational wave observatory, where the TM is enclosed 
by a rectangular capacitive sensor electrode housing - re- 
ferred to here as gravitational reference sensor (GRS) - 
with 3-4 mm gaps||. In the first configuration (1TM), 
the cubic TM is suspended on axis and thus a pure rota- 
tional damping is measured. In the second configuration 
(4TM), four TM are each displaced by roughly 11 cm 
from the pendulum rotation axis in a cross configuration. 
This results in both translational and rotational damp- 
ing inside electrode housings that surround two of the 
TM (one housing is the GRS mentioned above and the 
other is a larger housing with 6-8 mm gaps). Measure- 
ments with the 4TM were also performed after removing 
the GRS housing (labeled "w/o GRS" in the plot), in or- 
der to isolate its proximity effect on the total damping. 
From each dataset the zero-pressure background pendu- 
lum damping has been subtracted off using a linear fit. 
This is nearly irrelevant for the 1TM data, which employs 
a low damping fused silica fiber [l4[ - Q ~ 10 6 - while for 
the 4TM experiment, using tungsten fibers with Q « 
2000 - 3000, the fiber damping is ten times larger than 
the gas damping contribution at the lower pressures of 
the study. Here, uncertainty in the residual damping is 
the dominant source of error and causes the large spread 
in the data (several points with error bars overlapping 
zero have been omitted from the log-scale plot, but have 
been used in obtaining the fits to the 4TM data). 

The UTN data demonstrate the proportionality of gas 
damping to pressure, for more than three decades in pres- 
sure in the 1TM study and more than a decade in pres- 
sure in the 4TM study, with the wide dispersion men- 
tioned above limiting the validity at lower pressures. All 
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FIG. 6: Plot of experimental data from the UTN gas damp- 
ing experiment [5j featuring cubic test masses in two differ- 
ent torsion pendulums dominated by rotational (1TM) and 
translational (4TM) squeezing, the latter with two different 
conditions of proximity to the surrounding apparatus. The 
pressure independent damping has been removed from each 
dataset to show only the gas contribution. Linear fits to each 
dataset are shown. The simulation predictions for the three 
cases are shown as thick black lines, while the infinite volume 
limit predictions for the two pendulums are shown in dashed 
gray. 



the data demonstrate a clear - factor 3 to 6 - excess 
above the predictions of the infinite volume model calcu- 
lated from Ref . [J] . Additionally, the measurable decrease 
in gas damping for the identical pendulum upon removal 
of the GRS in the 4TM pendulum can only be attributed 
to proximity damping. The three datasets are consistent, 
at the 20% level, with simulations that employ the same 
technique as that described for the LIGO configuration 
in Sec. [TTl with the torque damping coefficient calcu- 
lated from the simulation torque fluctuations. Though 
based on the same basic cubic geometry and gap sizes, 
the different measurements probe both the translational 
and rotational aspects of squeezing, which have slightly 
different dependence on gap. 

Not shown here, the UTN experiments [5| also verified 
the Brownian nature of the force fluctuations associated 
with the gas damping, detecting a frequency independent 
increase in the pendulum torque noise floor consistent 
with S N = 4:k B Tf3 rot . 

The University of Washington (UW) experiments @ , 
with results shown in Fig. measure the gas pressure 
ringdown damping of a two different torsion pendulums. 
Both inertial members, labeled NTA and LTA, are thin 
rectangular plates in close proximity to a second parallel 
plate, with the NTA dimensions 42 mm x 31 mm x 2 mm 



FIG. 7: Plot of experimental data from the UW gas damp- 
ing experiment^ for two different torsion pendulums, la- 
beled NTA and LTA, both with a rectangular plate geometry 
in proximity of a second parallel plate at distance d. The 
pressure independent damping has been removed from each 
dataset to show only the gas contribution. Simulation points 
are shown with black x , with infinite volume limits shown in 
dashed gray. Dotted lines connecting simulation points are 
included as a guide for the eye. 



and the LTA dimensions 114 mm x 38 mm x 450 /im. 
Separations from 100 /im to 1 cm are studied, covering a 
range of aspect ratios - b/d, with b the shorter dimension 
of the rectangular plate - of 4 to 400, thus well brack- 
eting the LIGO geometry, where 2R/d w 70. The data 
have been normalized for linear pressure dependence, as 
well as for the different atomic masses of the residual gas 

1 /2 

species, mostly N2 and H2O, using p ~ pm . For each 
dataset the pressure independent damping has been es- 
timated, and removed, using the zero-pressure intercepts 
from linear fits to groups of damping measurements per- 
formed at fixed gap d. 

In addition to distributing statistically around a sin- 
gle curve following normalization for pressure, the data 
demonstrate the increase in damping with decreasing 
gap, by roughly two orders of magnitude for a factor 30 
decrease in gap for the LTA data. Also shown with the 
data in Fig. [7] are the predictions from the torque noise 
simulation performed for this geometry as in Sec. [TXJ 
The simulation /3 rot agrees with the experimental data to 
within 30% across two decades of d, with residual differ- 
ences considered compatible with systematic and statis- 
tical uncertainties in assigning pressures and gas species. 

Even for the largest gaps studied, b/d « 4 for LTA, 
excess proximity damping is dominant over the infinite 
volume contribution. The simulation for large d is com- 
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plicated by the presence of additional fixed control sur- 
faces, adding proximity damping that is independent of 
riflEl]. As seen in Fig. inclusion of geometrically simple 
control electrode surfaces in the LTA simulation create a 
relevant increase of the large d predictions above the infi- 
nite volume prediction [30( | . Thus, while the experimental 
data are consistent with the simulation at total damping 
values only several times larger than the infinite volume 
limit, they do not allow a quantitative comparison with 
the infinite volume predictions. 

The UW results are incompatible with an elastic, spec- 
ular scattering description of the molecule - wall inter- 
action. With elastic scattering, the molecule velocity 
parallel to the plate surfaces is unchanged as it bounces 
through the gap, and so the effective escape time r be- 
comes independent of gap. Following Eqn. [Til this would 
give force noise S ~ (3 ~ dT 1 . The steeper gap depen- 
dence observed in Fig. with logarithmic slopes of ap- 
proximately — 1.5 in the decade of smallest gap in both 
configurations, is thus incompatible with an elastic scat- 
tering scenario. Along with the success of the simulations 
in explaining the data, at least for the systems under 
study here - involving gold coated surfaces and mostly 
water vapor and nitrogen residual gas at room temper- 
ature - the diffuse inelastic scattering hypothesis is well 
verified. 

Taken as a whole, the two experimental campaigns 
measure proximity-enhanced gas damping at levels that 
agree with simulation predictions to better than 30% over 
several orders of magnitude in both pressure and gap size. 
Such a quantitative test of the physics and simulation 
techniques employed here for Advanced LIGO, in similar 
macroscopic systems, strengthens the confidence in the 
predictions for the gravitational wave sensitivity, which 
are discussed in the next section. 



V. IMPACT ON GRAVITATIONAL WAVE 
DETECTION CAPABILITIES AND DISCUSSION 

The data for Sp and r with d = 5 mm, in Figs. [3] 
and |4] correspond to the current Advanced LIGO design 
geometry. The low frequency noise power increase (or 
equivalently, the increase in the gas damping coefficient 
/3), with respect to the infinite volume limit is near 40, or 
slightly more than 6 in linear force noise spectral density. 

Considering a design gas pressure of 2xl0 -6 Pa (1.5 
x 10 -8 Torr) of H 2 , this contributes a white low frequency 
force noise S 1 /* 20 fN /Hz 1/2 , or an acceleration noise 

of nearly 0.5 fm/s 2 /Hz 1 ^ 2 . This acceleration noise is as- 
sumed to act equally on each of the four TM in the two 
interferometry arms. Converted into an equivalent gravi- 
tational wave strain noise with 4 km interferometry arms, 
the white force noise assumes an approximate /~ 2 shape 

1/2 

in S h ' at low frequencies, and, as shown in Fig. [51 rivals 
quantum fluctuations as the leading noise source in the 
10 - 30 Hz band [3,31]. This can have a significant im- 
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FIG. 8: Gravitational wave strain measurement noise contri- 
butions for Advanced LIGO Dashed and solid blue curves 
correspond to the noise contributions from gas damping noise 
with 2 x 10 -6 Pa of H2 for the nominal geometry (all 5 mm 
gaps between all four TM - reaction masses pairs) and the al- 
ternate geometry with a 20 mm gap adjacent to the two input 
TM. The total measurement noise in the two geometries are 
shown with thick gray curves. 

pact on the detection horizon for observing neutron star 
inspirals, even more so for the merger of more massive 
black hole sources, which have more signal contribution 
at the lowest frequencies. 

We note that the calculated diffusion time is fast 
enough, 700 [is, that the corresponding 3 dB rolloff fre- 
quency of roughly 230 Hz is above the band where gas 
damping is critical to the detector sensitivity, with pho- 
ton shot noise limiting sensitivity at these frequencies for 
the advanced detectors. 

The non-trivial impact of this noise source on the Ad- 
vanced LIGO sensitivity has raised the attention of the 
collaboration [It]], and we consider here a few options for 
mitigation. For the central input TM (ITM), increasing 
the gap between ITM and the reaction mass to 20 mm is 
feasible, which reduces the associated force noise power 
on the ITM by roughly a factor 7 (see Figs. [3] and [5J. 
The ITM reaction masses, called compensation plates, 
are used for correcting the beam wavefront for thermal 
lens effects and for relatively low force authority elec- 
trostatic actuators. The reaction masses associated with 
the distant end TM (ETM), on the other hand, are re- 
sponsible for higher authority electrostatic force control 
of the interferometer cavity length, and increasing the 
gap would create a serious limit on the maximum possi- 
ble force. As such, the current baseline design maintains 
5 mm gaps for the ETM. Mitigation for the ETM force 
noise could be achieved by improving the vacuum pres- 
sure at the end stations or, if necessary, adopting a more 
complicated, possibly ring-like, geometry for the electro- 
static actuation reaction masses. The gas damping and 
total gravitational wave strain with the increased ITM 
gaps is also shown in Fig. [5J 

It is worth reconsidering one initial assumption used 
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in this analysis, namely that the steady state pressure in 
the gap is equal to that of the surrounding gas reservoir. 
In practice, any local outgassing from the surfaces inside 
the gap will create a net radial molecular flux and thus 
a local pressure increase above that of the surroundings. 
This increased effective pressure in the gap increases the 
corresponding force noise, and thus maintaining the pre- 
dictions obtained for 2 /zPa of H2 limits the tolerable 
outgassing from the test and reaction mass surfaces. The 
pressure increase from a total gas outgassing flux Q (mea- 
sured in [Pam 3 /s]) can be estimated from the macro- 
scopic flow estimates of Sec. IIII A[ 

Ap = qz m = Q^r d • (22) 

Given the estimated values of r in Sec. [TTl maintaining 
Ap well below 2 /zPa in the 5 mm gap will require out- 
gassing well below 10 -6 Pa m 3 /s. Per surface area, this 
corresponds to 10 -8 mBarl/s/cm 2 , using more familiar 
vacuum technology units. This represents an important 
- but, with reasonable pumping times and temperatures, 
feasible - requirement on outgassing performance, which 
is likely to be dominated by water outgassing from the 
optical surfaces fl8l [l9j| . 

In closing, we discuss a series of design concerns rele- 
vant to other observatories with goals similar to those 
of Advanced LIGO, such as Advanced Virgo [2(| and 
the Lar ge-s cale Cryogenic Gravitational wave Telescope 
(LCGT [2"l[ I22I] ) . and to more ambitious future facilities, 
such as the Einstein Telescope (ET [23j]). The first issue 
is geometrical. Instead of simply increasing the distance 
between similar shaped test and reaction masses as in 
Advanced LIGO, which improves the linear noise spec- 
tral density roughly proportionally to d, the proximity 
enhancement can be effectively removed for a reaction 
mass in the form of an external coaxial ring around the 
cylindrical TM. In this geometry, TM motion along the 
cylinder (and GW-measurement) axis does not compress 
the gas between the TM and reaction ring. As such there 
is no gas flow in the gap and no related excess proximity 
dissipation. The same physics explains the absence of 
excess force noise perpendicular to the cylinder axis in 
the Advanced LIGO geometry seen in Sec. ILT1 

Advanced Virgo employs such an external ringlike re- 
action mass. There is an additional cylindrical thermal 
compensation plate facing the TM cylinder, but at a dis- 
tance of 20 cm [2(| and thus virtually irrelevant for its 
proximity (see Fig. [3]). Considering the vacuum require- 
ments, dominated by 10~ 7 Pa of H2O, and TM dimen- 
sions similar to those of Advanced LIGO, one expects a 
number several times below that given in Eqn. [3l and 
thus not a critical issue for the target sensitivity. 

Proximity concerns aside, gas damping can be reduced 
by acting on the factors relevant to the infinite volume 
limit, concerning the TM thermal vacuum environment, 
which enter as a factor p^^^T^-^m^ 4 in linear spectral 
density, and the TM geometry, entering roughly as the 
ratio of the square root of the the surface area (R) to mass 



(R 2 h). LCGT will employ cryogenically cooled TM, at 
20 K, with expected pressures near 2x 10~ 7 Pa dominated 
by H2. LCGT will thus enjoy a rough factor 6 noise im- 
provement over Advanced LIGO before considering ge- 
ometry (LCGT will also likely employ a ring geometry 
for the reaction masses (24j). Though the TM are slightly 
smaller than those of Advanced LIGO (30 kg of sapphire 
rather than 40 kg of fused silica), there should be a net 
improvement with LCGT and thus gas damping should 
not be critical for their design sensitivity. 

Several cryogenic issues complicate the gas damping 
analysis. The molecules striking the TM are assumed 
to have been thermalized by a cryogenic thermal shield, 
but some alteration of the velocity distribution by ther- 
mal gradients is still possible. Low temperatures should 
enable lower pressure with cryopumping, freezing out all 
relevant species except H 2 . More delicate, but potentially 
helpful, is the question of H2 sticking on the TM. Stick- 
ing times for H2 on cold solid surfaces rapidly lengthen 
with decreasing temperature, with a strong substrate de- 
pendence, from the [is - ms range at 20 K to the 1 s 
- 1 year range at 10 K (see App. |A"|) , While we do not 
have published absorption parameters for H2 on the most 
probable oxide TM coatings, such as Si02 and Ta205, it 
is likely that sticking times can rival or exceed time-of- 
flight to surrounding surfaces (of order 20 /is for 5 mm 
at 20 K). This lengthens the relevant escape time t and 
progressively lowers the maximum frequency for which 
proximity effects are relevant, even to below the band of 
interest [32j|. 

The Einstein Telescope, which has a dedicated low fre- 
quency interferometer, aims at sensitivities of roughly 
10- 24 /Hz 1/2 at 10 Hz with an increased armlength of 
10 km and thus requires a factor 30 improvement over 
the Advanced LIGO acceleration noise. This will require 
both an improved geometrical configuration, such as the 
just discussed ring reaction mass design, and some combi- 
nation of improved thermal vacuum conditions and larger 
TM. The specified conditions (10~ 8 Pa and 10 K) alone 
give a factor 30 improvement over the same values for 
Advanced LIGO. TM of roughly 200 kg and 500 mm di- 
ameter will give an extra factor 3 improvement for the 
infinite volume prefactor. Assuming that the target pres- 
sure is reached, gas damping does not appear to be a 
dominant noise source for ET. 



Appendix A: Estimation of molecular sticking time 

The sticking time, which we will call r s t, or average 
time that a molecule spends on the surface of the test 
mass or adjacent apparatus between impact and reemis- 
sion, can be evaluated via its influence on the two dimen- 
sional equation of state of the adsorbed film in equilib- 
rium with its gas vapor. This is characterized by number 
density n^D {p,T). In equilibrium, the rate of molecules 
(per area) that desorb from the surface to enter the gas 
phase, ^^j- , must be equal to the number of gas molecules 
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that strike the surface (25|. Integration over the Maxwell- 
Boltzmann velocity distribution gives the latter, and we 
findfH 



(2nm k B T) 



1/2 



n 2 p 

Tat 



(Al) 



For the low gas pressures and relatively high tempera- 
tures of interest to our experimental system, we consider 
the adsorbed molecules as an ideal non-interacting 2D 
gas in the xy plane, with a substrate potential V (z) in 
the dimension perpendicular to the surface. Approximat- 
ing the surface potential with a well depth D and elastic 
constant toqw 2 , this gives energy levels 



h 2 k 2 
2mo 



where k represents the 2D wave vector. In the classical 
limit of low occupation number, this yields the following 
equation of state 
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Combining this with Eqn. IA11 we obtain 
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An intuitive, semi-classical interpretation of this sticking 
time emerges in the limit of classically excited oscillator 
states, with Tllo <C fcgT, which yields 
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This limit corresponds to a molecule oscillating around 
the bottom of the potential well, "attempting escape" 
once per classical oscillation, with an Arrhenius escape 
probability of exp ^4yr . 

Substituting a harmonic well for the true Van der 
Waals potential with Coulomb repulsion is an approxima- 
tion, both for the level spacing and for the finite number 
of bound states. Additionally, this model treats possible 
rotational excitations of the molecule as equally popu- 
lated in the gas and adsorbed phases. As such, we use 
this model only to have an order of magnitude estimate 
of sticking times, to give a comparison with the typical 
times of flight between collisions on adjacent test mass 
and reaction mass surfaces. 

Typical measured values for the adsorption well depth 
of H2 on various solid substrates, including both di- 
electrics and metals, range from 300 to 600 K, with os- 
cillator level spacings of order several hundred K [26| . 
This results in room temperature sticking times of order 
1 ps. The distinction between the classical and quantum 
excitation levels of harmonic oscillator is a minor correc- 
tion to the sticking time estimates here, and the sticking 



times correspond to several to tens of classical oscillation 
cycles. More easily polarized molecules such as N 2 , CH 4 , 
and H2O can find adsorption well depths of order 1000 K 
with slightly lower vibrational energies, resulting in stick- 
ing times of order 10-100 ps. For all realistic residual gas 
species, the room temperature sticking times are sub-ns, 
well below the typical times for a gas molecule to cross 
a 5 mm gap, which is of order 5 /is for H2 and longer 
for larger molecules. As such, the sticking time of resid- 
ual gas molecules in Advanced LIGO is not expected to 
create a relevant departure from the time dependence of 
force fluctuations given using the approximation of im- 
mediate reemission of incoming gas molecules. 

For cryogenic temperatures, such as the 20 K fore- 
seen for LCGT, the situation is quite different. Bound 
molecules are frozen into the ground state of the adsorp- 
tion potential well, with a binding energy (D — 4^) . The 
sticking time for H2 on the dielectric surface of foreseen 
oxide coatings such as SiC-2 or Ta20s becomes very sen- 
sitive to the adsorption parameters - for which we do not 
have published values - and the exact temperature. For 
instance, for H2 on the oxide MgO, with binding energy 
330 K, the estimated sticking time is of order 0.5 ms, 
with a 10% change in either binding energy or temper- 
ature producing a factor 5 change in the sticking time. 
For H2 on NaCl (280 K binding energy) and on graphite 
(440 K), the sticking times are 3 /is and 7 ms. At 10 K, 
the range of T s t for these species becomes 10 s to 1 year. 

In the Advanced LIGO-like geometry at cryogenic tem- 
peratures, it is thus conceivable that the sticking time 
will be longer than the time for molecules to cross the 
intercylinder gap. Though this complicates the analy- 
sis, the general effect on gas damping can be understood 
with the force shot noise model discussed in Sec. IIII CI 
At a given p and T, sticking will not change the total 
momentum exchange associated with a molecule's stay 
inside the gap or the rate at which molecules enter the 
gap, but it will lengthen the typical time spent inside the 
gap by N T s t (N is the typical number of collisions that 
a molecule makes before diffusion out of the gap, Sec. 
IIII1 of order 100 for 5 mm gaps in Advanced LIGO). As 
such, the low frequency force noise should be unchanged 
by sticking, but the longer diffusion or correlation time 
results in a lower frequency cutoff for the excess noise 
from any proximity effects. This could push any excess 
proximity damping to frequencies below the band of in- 
terest. 



Appendix B: Calculation of (Aq 2 ) 

For a zero mean, statistically stationary random vari- 
able F (t) with mean defined over some time interval 
[0,T ], F = tjt J Q F (t) dt, the expected fluctuations in 
F can be calculated in terms of the autocorrelation func- 
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tion Cff (5) = (F (t) F{t + 6)), 



(F 2 ) = jr 
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In our case, F (t) is the instantaneous force felt by the 
TM due to molecular impacts. 

The single-sided power spectrum Sf (w) is defined 
by the Fourier transform of Cff (S), Sf (w) = 
2 J^ 00 Cff (S) exp (— iu8) dS. For a process character- 
ized by white noise passed through an effective low-pass 
filter, such as our approximation of the gas damping sys- 
tem studied here, with Sf = ; , r ,i , the corresponding 

1 + (cjt) 7 ^ ° 

autocorrelation function can be calculated to be 
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FIG. 9: Illustration of the geometry used for calculating the 
mean square lateral displacement and the mean time of flight 
for a molecule emitting into the gap from the center of one 
of the cylinders. The angle 8c represents the critical value 
of the polar angle 9, beyond which emission from the center 
results in the molecule escaping the gap laterally in a single 
flight. 



Inserting this result into Eqn. IB1I and integrating, we 
obtain 
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As the total momentum exchange with the TM in a mea- 
surement time T is simply Aq = FT , we recover the 
result in Eqn. [4] 



Appendix C: Calculation of (r 2 ) and (i) 

Estimates for (r 2 ) and (t) are obtained by direct in- 
tegration over the distribution of outgoing molecular ve- 
locities for an emitted molecule 
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where 9 and cj> are the polar and azimuthal emission an- 
gles. For both calculations, we divide the distribution 
into a range of polar angles 9 (see Fig. for which a 
molecule emitted from the center of a TM face strikes 
the opposing cylinder (9 < 9 C ) and another for which 
the molecule escapes the gap directly (9 > 9 C ), with the 
critical angle defined by the aspect ratio, tan6* c = -j- 
For the average time of flight, we obtain 
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For the calculation of (r 2 ) only the angular distribution 
P (9, <h) dil = — cos 9 dfl is relevant. We note that those 



molecules, emitted for 9 > 9c-, that directly leave the 
gap on a single jump from the center are weighted with 
a lateral displacement that is "saturated" to be R. We 
obtain 
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We note that while the average time of flight remains 

1 /2 

finite for cylinders of infinite radius - (t) — > ) ^7, 
as we use in the approximate calculation (Eqn. 1121) 
the mean square displacement diverges, albeit slowly, as 
R/d — > 00. Saturating the maximum displacement to R 
makes physical sense and corresponds to the case of es- 
cape on a single bounce. The logarithmic dependence of 
(r 2 ) is a fundamental difference with the result of Ref. 
[8|, where the relevant random walk step size is taken 
to be (r) w ^ (and not the RMS ^/(r 2 ) as used here). 
This would result in analogous expressions for N, t, and 
Sf similar to those expressed here (Eqns. [TBIfTSj) but 
without the logarithmic factor. The logarithmic factor, 
also relevant to similar expressions for the molecular flow 
impedances of short pipes (for example Ref. 27]) is re- 
sponsible for the nearly power law dependence of j3 (or 
Sf) on 1/d, with an exponent that is not exactly 2, but 
rather approaches 2, from below, only for vanishing gap, 
as observed here and in the three other numerical and 
experimental studies cited here (Refs. H, H, §). 
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